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Abstract 

We review the theory of e'/e and present an updated phenomeno- 
logical analysis using hadronic matrix elements from lattice QCD. 
The present status of the computation of e'/e, considering various 
approaches to the matrix-element evaluation, is critically discussed. 



* based on the talks given by M.C. at "Les Rencontres de Physique de la Vallee d'Aoste", 
La Thuile (Italy), 27 February-4 March 2000 and by G.M. at the "XXXVth Rencontres 
de Moriond", Les Arcs 1800 (France), 11-18 March 2000. 



1 Introduction 



The latest-generation experiments, aiming to obtain s'/e with a 10 4 accu- 
racy, measured up to now 

e^_{ (28.0 ± 4.1) x 10" 4 KTeV [1] 

e ~ j (14.0 ±4.3) x 10- 4 NA48 [| ^ ' 

By combining these results with previous measurements, the latest world 
average reads 

- = (19.3 ±2.4) x 10~ 4 , (2) 

which is definitely in the 10~ 3 range. Given the differences in the results of 
eq. (]]]), the quoted error is, however, debatable [[|. 

On the other hand, theoretical estimates in the Standard-Model typically 
correspond to central values in the 10 -4 range although, given the large 
theoretical uncertainties, values of the order of 10 -3 are not excluded. The 
explanation of the difference between SM predictions and experimental values 
calls either for some missing dynamical effect in the hadronic parameters or 
for physics beyond the Standard Model. In the last few months, several 
studies exploring both possibilities have been published. 

In this paper, the theoretical status of e' / e is reviewed and updated results 
obtained by using (whenever possible) hadronic matrix elements computed 
with lattice QCD are presented. Other theoretical approaches and recent 
attempts to "improve" the accuracy in the determination of the hadronic 
matrix elements (mostly to improve the agreement between theoretical esti- 
mates and measurements) are also discussed. 



2 Basic formulae 

Direct CP violation, occurring in K° decays, is parametrized by e'. In terms 
of weak-Hamiltonian matrix elements, this quantity is defined as 
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are the eigenstates of the CPT-conserving Hamiltonian describing the K°-K° 
system, namely 

rr M V ( M Mi 2 \ « / r r 12 \ 

We introduce the isospin amplitudes 

Aje iSl = (nn(I)\H w \K°) , A*e iSl = (nn(I)\H w \K°) , (6) 

where, in virtue of Watson's theorem, the 5js are the strong-interactions 
phase shifts of tttt scattering. In the approximation ImAo -C KeAo, Im^ <C 
ReA 2 and to = Rev4 2 /Rev4o <C 1 (the latter coming from the A J = 1/2 
enhancement in kaon decays), one finds 

.ImAo _, e in/4 ( ImMi2 + ImA c 
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e' ~ i- — =— Im — ~ i =— - - (a; ^irn^ - Iim4 ) . (7) 
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Using the experimental value 0] 

Al ge' = | + 5 2 - 5 = (48 ± 4)° « tt/4 , (8) 

one finally gets 



where the last expression includes isospin breaking contributions due to tt-i] 
mixing encoded in Q IB (A' 2 = A 2 — ujVL ib Aq) ||. In the prediction of e'/e, 



2 



uj and ReA are taken from experiments, whereas lmA 2 are the computed 
quantities. 

The calculation of the real part of the amplitudes, and hence of a;, is 
one of the longest-standing problems in particle physics: in spite of several 
decades of efforts, nobody succeeded so far to explain the AI =1/2 rule in 
a convincing and quantitative way. The calculation of lmA and Iul4 2 is 
of comparable difficulty Since the imaginary parts entering e'/e, however, 
are not directly related to the real ones, and the operators of the effective 
Hamiltonian contribute with different weights in the two cases, it is conceiv- 
able that Iul4 and Im^ be computed in spite of the difficulties encountered 
in calculations of the AI =1/2 rule. On the other hand, as discussed below, 
one cannot exclude some common dynamical enhancement mechanism which 
produces large values of both ReA and e'/e. 



3 e'/e in the Standard Model 

The natural theoretical framework in dealing with weak hadronic decays is 
provided by the effective Hamiltonian formalism. Indeed the operator prod- 
uct expansion allows the separation of short- and long-distance scales and 
reduce the problem to the computation of Wilson coefficients, performed 
in perturbation theory, and to the calculation, with non-perturbative tech- 
niques, of /oca/-operator matrix elements. 

At the next-to-leading order (NLO) in the renormalization-group im- 
proved expansion, the 4-active-flavour (ra^ > /x > m c ) AS = 1 effective 
Hamiltonian, relevant for e'/e, can be written as 

C 1 (/i)(Q 1 (/i) - Q C M) + C 2 (/i) (Q 2 (jjl) - Q C M) 

+r^C i (/i)Q i (/i)} , (10) 

i=i 

where Gf is the Fermi constant, X q = V q dV* s and r = —\t/\ u (Yq l q ] being the 
CKM matrix elements). The CP-conserving and CP- violating contributions 
are easily separated, the latter being proportional to r. 



H W = (1-r) 
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Neglecting electro- and chromo-magnetic dipole transitions, the operator 
basis includes eleven independent local four-fermion operators. They are 
given by 

Ql = {s a d a )v~A{Uj3Up) V -A , Q2 = (Sad(3)v-A(U/3U a )v-A, 

Ql = {s a d a ) v _ A (cpCp)v-A Q2 = {s a dp) V -A{cpCa)v-A 1 

Q3,5 = (s a d a ) V -AY,q(<ll3qi3)vTA, Qifi = V-A Eg W/3<?a) V=Rl , 

(11) 

where (q a q'g)v±A — 9a 7^(1 ^ Ts)?^ a an d ft are colour indices, and the sum 
index q runs over {d,u, s,c}. The operator Q2 appears in the Fermi Hamil- 
tonian at tree level. The operators QztQq are generated by the insertion of 
Q2 into the strong penguin diagrams, whereas Q7-Q9 come from the electro- 
magnetic penguin diagrams. Both classes of operators are relevant for e'/e. 
Further details on the NLO AS = 1 effective Hamiltonian can be found in 
refs. f§. 

Using eq. fllCf ), one can readily express Aq and A2 in terms of matrix 
elements of the operators in eq. ([0]). It is customary to write 




(12) 

W 2 = -G F lm(v; s V td ){C 7 B 3 7 /2 (j-j)+C 8 Bl /2 (Y~j) 

In the previous equation, the relevant matrix elements are given in terms of 
i?-parameters defined as 

(7T7r(0)|Q,|K> = B^imrMlQiWviA, 

(7r7r(2)|Q i |i^> = B^innWlQ^vjA, (13) 
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where the subscript VIA means that the matrix elements are calculated 
in the vacuum insertion approximation. VIA matrix elements are given in 
terms of the three quantities 



Contrary to X and Z, Y does not vanish in the chiral limit, as a consequence 
of the different chiral properties of the operators Q7 and Q$. Moreover, 
whereas X is expressed in terms of measurable quantities, both Z and Y 
depend on the quark masses which must be taken from theoretical estimates. 

Note that some VIA matrix elements seems to show a quadratic depen- 
dence on the strange quark mass m s through Y and Z. This is true as 
long as one fixes the kaon mass to its experimental value and neglect the 
m s dependence of the B parameters. The apparent quadratic dependence of 
the matrix elements on m s has been exploited in ref. to claim that large 
values of e'/e can be obtained with suitably "small" strange quark masses. 
The actual dependence of the full matrix elements on m s is, however, very 
different. Indeed the ratio M^/m s is essentially independent of m s (up to 
small chiral-symmetry-breaking terms) since it corresponds to the value of 
quark condensate. This is explicitly verified in lattice calculations, where a 
strong correlation between the value of the strange quark mass used in the 
VIA matrix elements and the value of the corresponding B parameters is 
observed, so that the m s dependence in the physical matrix element almost 
cancels out ||. This is why one should always use 5-parameters and m s 
consistently computed together (e.g. in the same simulation on the lattice) 
or, even better, matrix elements given in physical units without any reference 
to quark masses || || . 



X 



U (Ml - Ml) 



Y 
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4 Hadronic matrix elements calculation 



Any prediction of e' / e must undergo the non-perturbative calculation of the 
relevant hadronic matrix elements. Theoretically, this calculations has to 
meet two requirements 

1. to be applicable up to perturbative energy scales; 

2. to keep under control the definition of the renormalized operators and 
their consistent matching to the Wilson coefficients. 

Failure to meet these requirements indicates that the method cannot achieve 
the necessary NLO accuracy, as often the case with phenomenological mod- 
els. Presently, however (or may be for this reason), experimental data are 
more easily accommodated by such models than by methods based on first 
principles, as lattice QCD. 

With this caveat in mind, we list and briefly comment on various ap- 
proaches that have been used in the literature. Various predictions of e'/e, 
obtained using different methods, are shown in the compilation of fig. [I], 
together with the present experimental world average. 

Lattice QCD [0 

In principle, Lattice QCD is the non-perturbative method for computing ma- 
trix elements. Being a regularized version of the fundamental theory, it allows 
a complete control over the definition of the renormalized operators, both at 
the perturbative and non-perturbative level. In addition, present simulations 
use inverse lattice spacings of 2 GeV or larger and therefore the perturbative 
matching with the Wilson coefficients can be safely performed. Indeed, by 
using non-perturbative renormalization techniques, the matching scale in lat- 
tice calculations could be pushed to values as large as ~ 10 GeV. Although, 
so far, these methods have only been implemented in the calculation of the 
strong coupling constant 0] and of the quark masses [O], they will certainly 



be extended to the four-fermion operators of the weak effective Hamiltonian. 
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Figure 1: Compilation of recent theoretical predictions for e'/s. 



When this will be the case, the error in the matching procedure will become 
negligible. 

For many years, a general no-go theorem |TB|] of Euclidean field theory 
has prevented the direct extraction, in numerical simulations, of the physical 
matrix elements with more than one particle in the final state. For this rea- 
son, present lattice determinations of (7nr\Qi\K) are obtained from (n\Qi\K) 
using lowest-order chiral relations (i.e. soft-pion theorems). This means that 
final-state interactions are not taken into account and that large chiral cor- 
rections may be present |14| Q In addition, only some of the matrix elements 
needed for computing e' / e are presently available. In particular, the Qq ma- 
trix element, which is expected to give the most important contribution to 
e'/e, has not been successfully computed yet. 

Several theoretical progresses have opened a window of opportunity this 

1 Other typical lattice systematics, such as the quenching, are expected to play a lesser 
role in this context. 
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year. In the past, a proposal to circumvent the no-go theorem of ref. | T3[ 
was made. The main idea was to extract the relevant matrix elements by 
studying suitable Euclidean Green functions at small time distances |T5| . The 
weakness of this method, however, was that it relied on mo del- dependent 
smoothness assumptions which could lead to uncontrolled systematic errors. 



A big progress has been made in ref. ||16| , where it was rigorously proven 
how to relate the matrix elements extracted on a finite volume in lattice 
simulations to the physical (7nr\Qi\K) amplitudes. Moreover, it has been 
shown that the smoothness hypothesis of ref. |~5| is unnecessary and that 



the physical K — > nn matrix elements can be, at least in principle, extracted 
from Euclidean correlation functions at finite-time distances |17|| . Although it 



will take some time before these approaches will be implemented in practice, 
they certainly open new perspectives to lattice calculations. 

More details on the present status of lattice matrix elements, and pos- 
sible developments in the near future, will be given in the discussion of the 
phenomenological analysis in the next section. 

Phenomenological Approach ||18|| 

In this approach, one basically attempts to extract information on the ma- 
trix elements relevant for e'/e by combining the measured values of the CP- 
conserving amplitudes with relations among different operators that can be 
established below the charm threshold under very mild assumptions (for de- 



tails see ref. [fL9|1 ). This procedure can be performed consistently at the NLO, 
allowing the extraction of matrix elements of well-defined renormalized op- 
erators. 

Unfortunately, the leading contributions to e'/e, namely the matrix el- 
ements of Q$ and Q 8 , cannot be fixed in this approach. Moreover, the 
method only works below the charm threshold where higher-order perturba- 
tive and power corrections (in l/m c ) may be large. In practice, for these 
matrix elements, Buras and collaborators have always used inputs coming 
from other theoretical sources, in particular lowest-order 1/N expansion or 
lattice calculations. 
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Chiral+l/iV Expansion Q21J] 

This method relies on the non-perturbative technique originally proposed by 
Bardeen, Buras and Gerard |2l|]. In principle, the approach can be derived 



from QCD and allows the computation of all the matrix elements needed 
for calculating e'/e in a consistent theoretical scheme. In this framework, 
the Dortmund group computed the relevant matrix elements including the 
subleading corrections in both the chiral and the 1/N expansion 0. 

This approach suffers, however, from the presence of quadratic diver- 
gences in the cutoff that must be introduced, beyond the leading order, in the 
effective chiral Lagrangian. The quadratic cutoff dependence, which appears 
in non-factorizable contributions, makes it impossible a consistent matching 
between the operator matrix elements and the corresponding Wilson coeffi- 
cients, which depend only logarithmically on the cutoff. One may argue that 
the quadratic divergences will be cured and replaced by some hadronic scale 
in the full theory, which includes excitations heavier than the pseudoscalar 
mesons. In practice, since it is impossible to include the effects of higher- 
mass hadronic states, the cutoff is replaced with a scale of the order of 1 GeV, 
which is an arbitrary, although reasonable, choice. Since the divergent terms 
gives very large contributions to the matrix elements entering AI = 1/2 
transitions and e'/e, this introduces an uncontrolled numerical uncertainty 
in the final results. 



Chiral Quark Model [£2 



The xQM can be derived in the framework of the extended Nambu-Jona- 
Lasinio model of chiral symmetry breaking [p^j . It contains an effective 
interaction between the u, d, s quarks and the pseudo-scalar meson octet 
with three free parameters, two of which can be fixed using CP-conserving 
amplitudes. The Trieste group computed the 0{p A ) corrections to the rel- 
evant operators and found a correlation between the CP-conserving and 
CP-violating amplitudes so that, once the parameters of the model are fixed 



2 Indeed, some of the higher-order terms have been computed in the chiral limit only. 
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to provide the required octet enhancement, it is possible to predict e'/e. The 
nice feature of this model is that, to some extent, it accounts for higher-order 
chiral effects, which are not easily included, for instance, in lattice calcula- 
tions. The disadvantage is that the model dependence of the results can 
hardly be evaluated or corrected. 

Theoretically, this approach shares some of the problems of the 1/N ex- 
pansion mainly those related to the presence of quadratic divergences. These 
do not appear explicitly in the calculations of the Trieste group simply be- 
cause dimensional regularization is used. It remains true, however, that the 
scale and scheme dependence of the renormalized operators is not under con- 
trol at NLO. In order to deal with this problem, a third parameter of the 
model is fixed by imposing a sort of numerical '^"-independence to the phys- 
ical amplitudes. This recipe, while suggesting that some degree of "effective" 
renormalization-scheme independence can be achieved, has unfortunately no 
sound theoretical basis. Finally the correlation between the AI = 1/2 am- 
plitude and e'/e is subject to potentially large uncertainties for the following 
reason. The parameters necessary to estimate the matrix element of Qq are 
fixed by fitting the A J =1/2 amplitude. For this quantity, the contribution 
of Qq is rather marginal, whereas Qi and Q2 dominate. Thus any small 
uncertainty in the dominant terms, due for example to unknown 0(p 6 ) cor- 
rections ((9(p 4 ) corrections to the AI = 1/2 amplitude are of (9(100%)), may 
change drastically the determination of the matrix element of Qq, which is 
the dominant term for e' / e. 



Extended Nambu-Jona-Lasinio Model PB| 



An extended Nambu-Jona-Lasino model has been used in ref. |23| to com- 
pute the AI = 1/2 K — > 7r7r matrix elements and e'/e. The remarkable 
feature of this computation is the high order in the momentum expansion 
reached by the Dubna group. All matrix elements have been computed to 
0(p 6 ) and a good stability of the results has been found. In this respect, 
this approach is safer than the xQM. However it shares with the xQM all 
the other theoretical flaws mentioned above, and particularly the problem of 
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matching the short-distance calculations, since it is unclear which renormal- 
ized operators the amplitudes computed with the Dubna superpropagator 
regularization method correspond to. 



Generalized Factorization [| 

Generalized factorization has been introduced in the framework of non-leptonic 
B decays in order to parametrize the hadronic matrix elements without a- 
priori assumptions ||25|| . The basic idea is to extract from experimental data 
as much information as possible on the non-factorizable parameters. When 
needed, the number of independent parameters can be reduced using flavour 



symmetries, dynamical assumptions, etc. In ref. [23] the procedure has been 
applied to K — > tctc matrix elements. Unfortunately, in this case, the num- 
ber of independent channels that one can use to fix the parameters is small 
(essentially only the two CP-conserving amplitudes). For his predictions, 
the author of ref. [^3J was forced then to reduce the number of parame- 
ters by several "simplifying" assumptions, which are, however, questionable. 
Many parameters related to different operator matrix elements and to dif- 
ferent colour structures were arbitrarily assumed to be equal. In such a way, 
a correlation between CP-conserving and CP- violating amplitudes was ob- 
tained, but the final results crucially depends on the assumptions, which are 
hardly justifiable theoretically, and cannot be tested phenomeno logically in 
processes different from e'/e. 

a Models [|7|, |§ 

A possible mechanism to enhance the AI = 1/2 amplitude is the exchange 
of a scalar 1 = meson [27|. It also leads to an enhancement of e'/e, as 



recently studied in the framework of the linear [[7] and non-linear a 
models. While unable to achieve NLO accuracy, these models can produce 
the required correlation between the A J = 1/2 rule and e'/e, at least for some 
choice of the free parameters, such as the a mass. Also in this case, however, 
it is not easy to estimate the uncertainties and the model dependence of the 
theoretical predictions. 
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Other theoretical developments 

The marginal agreement between the SM predictions and the measured value 
of e'/e stimulated various attempts to "improve" the determination of the 
operator matrix elements, by including effects that were not considered pre- 
viously. In particular, new studies have been devoted to the calculation of 
isospin-breaking and final-state interaction effects. 

In most of the approaches isospin-breaking corrections are not included 
because they are beyond reach for these methods. These effects can be 
evaluated, however, a posteriori and included in the predictions. The leading 
effect in the chiral expansion is expected to come from 77-77-77' mixing, which 
can be computed following ref. ||. The resulting isospin-breaking effect is 
accounted for by the parameter Qjb which appears in eq. (|9|). Recently the 
calculation of Qjb has been updated by including the effect of ir° — 7] mixing 
at 0(p 4 ) P5[ . In addition, it has been pointed out in ref. |F2U| that new 
sources of isospin breaking appear, beyond the leading order, in the chiral 
Lagrangian. These terms may give large corrections to Qib- Unfortunately 
the calculation of the corrections is strongly model dependent and can only 
been taken as a warning on the potential importance of these effects. 

The problem of including final-state interactions is particularly relevant 
for lattice or lowest-order 1/N calculations, where rescattering effects are 
missing. It has recently been suggested that these effects could be included 
by using the measured tot phase shifts and dispersive techniques |L4 |. The 



resulting Aq amplitude would be enhanced by the inclusion of final state 
interactions, giving for e'/e a prediction much closer to the experimental 
value. These approach have been subject to several criticisms On 
the one hand the analytic structure of the considered amplitudes is unclear 
and the corresponding dispersion relations questionable. On the other, the 
computation of the dispersive correction factors, as derived in ref. [|T4 



is 



plagued by an irreducible ambiguity of the same order of the dispersive factors 
themselves. This uncertainty depends on the choice of the initial conditions 
which, as shown in JHJ, were arbitrarily chosen. For this reason, whereas 
final-state interactions are likely to give qualitatively a certain enhancement, 
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as argued in ref. the quantitative estimate of these effect is subject to 
very large uncertainties. As discussed in [P0|1 , lattice calculations could help 
in this respect by fixing the initial conditions in a unambiguous way. 

Some short, provocative comment is necessary at this point. If one could 
implement in the same calculation all the corrections which have been sug- 
gested to improve the accuracy in the determination of the matrix elements 
(low strange-quark mass, isospin-breaking effects, final state interactions, 
etc.), one would probably end up with a prediction of e'/e much larger than 
its experimental value. It is also quite astonishing that the effects which were 
not considered before, or those which have been revised in recent studies, all 
increase the theoretical prediction for this quantity and no one goes in the 
opposite direction. Finally, if the AI =1/2 rule and the large value of e'/e 
are a consequence of many effects which are all necessary in a conspiracy 
to give a large enhancement, it seems very unlikely that any of the existing 
theoretical approaches (including the lattice one), will ever be able to take 
them into account simultaneously at the necessary level of accuracy. 

During the completion of this paper, several new calculations of e'/e ap- 
peared: i) a new estimate of the Qq and Q^ 2 matrix elements using QCD 
sum rules has been presented jnj, with results for e'/e close to the exper- 
imental average; ii) within big uncertainties, very large values of e'/e have 
also been obtained using the 1/N and chiral expansions within the context 
of the extended Nambu-Jona-Lasinio model. In ref. [32], the proposal for 
controlling the scale and scheme dependence of renormalized operators using 
an intermediate X-boson has been implemented in the calculation of e'/e. 
We refer the reader to the original publications for more details. 



5 Results for e'/e using Lattice QCD 

In order to compute e'/e, besides hadronic matrix elements one needs the 
value of the relevant combination of CKM matrix elements, namely lmV t *Vtd- 
This is constrained by using the experimental information on |V^,|, \V u b\, 
Bd, s -Bd,s mixings and e combined with lattice results. Nowadays this has 
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become a standard way of determining the CKM-matrix parameters within 



the Standard Model, described for instance in refs. [[Hj, |33|. It is worth 
noting that the linear dependence of e'/e on ImV t *V t d is strongly reduced by 
the constraint on the CKM parameters enforced by the measured value of 
e. In the analysis reported in this paper, the same input parameters as in 



ref. IP- 0|| , with the exception of Qjb = 0.16±0.03 which is now taken from [28 
have been used. 



The discussion in ref. [JlOj about the current status of the lattice com- 
putation of the main matrix elements (7T7r\Qe\K) and (tttc\Q^ 2 \K) can be 
summarized as follows: 

• At present, the matrix element (tctt\Qq\K) is not reliably known from 
lattice QCD. The results with staggered fermions are plagued by huge 
corrections appearing in the operator renormalized using lattice pertur- 



bation theory p4| . Other attempts using Wilson fermions or domain- 
wall fermions were unsuccessful so far. 

The matrix element (7nr\Q$ \K) has been computed by several groups 
using different formulations of the lattice fermion action and different 
lattice spacing. A substantial agreement of the different determinations 
was found within 20% uncertainty. We use the value 

(HQ8 /2 (2GeV,MS-HV)|fO = 0.62 ± 0.12 GeV 3 (15) 

for the operator renormalized at fi = 2 GeV in the 't Hooft-Veltman MS 
scheme. This value corresponds to the B parameter B^ 2 = 0.93 ±0.18 
for a "conventional" quark mass m^ IS + m^ IS = 130 MeV at fi = 2 GeV, 



see eq. (|13| ) 



Using some reasonable assumptions for the less important contributions due 
to other operators, given that the largest uncertainty stems from our igno- 
rance of (inr\Q 6 \K) , a useful way of presenting the results is given by 



"(-25.3l|jl) GeV" 3 (nn\Q 6 (2 GeV, MS-HV) | K) - 6.3±J; C 7 



5 



x 10" 4 , 

(16) 
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where the matrix element of Q§ is considered as a free parameter. Notice 
that the the two terms in this equation are correlated and should not be 
varied independently. 

In order to compare eq. (|16|) with experiments, we have to make some 
assumption on the value of {txti\Qq\K) . We take the central value suggested 
by the VIA (or equivalently by the lowest-order 1/N expansion), namely 
B 6 = 1, with an uncertainty of 100%. This introduce a renormalization- 
scheme ambiguity, since VIA does not allow a proper definition of the renor- 
malized operators. For this reason, results obtained by taking two different 
central values for (itit\Q§\K) (corresponding either to B™(2 GeV)=l or to 
Bq DR (2 GeV)=l are presented, namely 

(HQ 6 (2GeV,MS-HV)|iQ = -0-4±0.4GeV 3 (Bf v = l±l), 
(7T7r|g 6 (2GeV,MS-HV)|ir) = -0.6 ± 0.6 GeV 3 (B^ DK = 1±1). ^ ' 
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In the two cases, we obtain 



'■i 

e'/e = (4.6±££ ± 0.4) x 10" 4 Monte Carlo 

-13 x 1(T 4 < e'/e < 26 x 10~ 4 scanning 

e'/e = (8.1^° 5 3 ± 0.3) x 10" 4 Monte Carlo ^ 18 ^ 

-13 x 10~ 4 < e'/e < 37 x 10~ 4 scanning 



The difference of the two results, contrary to what is often stated in the 
literature, is not the uncertainty associated to the renormalization-scheme 
dependence, but to a different choice of the value of the matrix element in a 
given scheme (the HV scheme in the example of eq. flT7])). At the NLO, the 
scheme dependence comes from higher-order corrections only and its effect is 
estimated by the second error given in eq. (0). The two figures of eq. ( |TB| ) 
correspond really to two different choices of the unknown value of (irn | Q e | K ) 
at a given scale (fj, = 2 GeV) and in a well-defined scheme (MS-HV). On 
the contrary, the two distributions of values for e'/e in fig. |2| include, for the 
same choice of (tttc\Qq\K) , two different ways to match Wilson coefficients 
and matrix elements for estimating the real scheme dependence due to higher- 
order terms. Both distributions refer to the case B 6 (2 GeV, MS-HV) = 1 ± 
1. The large error on the matrix elements of Qq obviously dominates the 
final uncertainty on e'/e and flattens these distributions. In ref. |18j a more 
optimistic error for the B parameter (B$ (2 GeV, MS-HV) = 1 ± 0.2) was 
assumed. 

6 Comparison with data 

Many of the Standard Model predictions shown in fig. [I] are below the present 
experimental world average. What does this imply? There are three legiti- 
mate answers: 

1. There is nothing wrong! For some specific choice of the input parame- 
ters, all the different approaches are able to reproduce to some extent 
the experimental data. In some cases the agreement seems to arise 
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naturally from the calculation PTl In other cases, this requires the 



adjustment of a few parameters [22] or a wise choice of several of them 



(often at the edge of the allowed range of values) [T3|, [TTJ. It is puzzling 
that most of the approaches, which suffer from intrinsic and irreducible 
uncertainties coming from the model dependence of the results, are in 



good agreement with the data. In the case of refs. [10|, [L8|, instead, 
this requires that all the quantities on which we have a poor control 
conspire in the direction to increase the theoretical value of e' / e. Thus, 
although unlikely in our opinion, the possibility that there is nothing 
wrong is not excluded. It may also well be that some of the models are 
indeed able to describe the underlying strong dynamics. 

2. There is something missing in the computation of the matrix elements. 
The long-standing problem of explaining the AI = 1/2 rule suggests 
that some important dynamical effect is at work in K — > tttt I = 
decays. Unfortunately, contrary to some old claim, there is no simple 
relation between the CP-conserving and CP- violating decays, which 
could explain the large value of e'/e on the basis of the enhancement 
of the AI = 1/2 amplitude. Indeed, it would be very interesting 
if a common dynamical mechanism could explain both of them. In 
terms of Wick contractions in the matrix elements, such a mechanism 
could be possibly provided by a large contribution from eye diagrams 
(aka penguin contractions) [|35| . From the lattice estimates, by taking 
B 6 as a free parameter, we can reproduce the experimental e'/e with 
P<? v (2GeV) ~ 2.4. As we have seen, all non-perturbative methods 
are affected by theoretical and/or computational problems which limit 
their accuracy. Among them, the models based on the chiral expansion 
also support the existence of some correlation between the AI =1/2 
rule and e'/e which is at least in qualitative agreement with the ob- 
servations. A possible exception is that of ref. p3 |. Thus we conclude 



that a real quantitative explanation is still to come. 
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3. Hadronic matrix elements are fine. New physics is at work. If the the- 
oretical calculations which gives low values for e'/e are correct, there is 
room for new physics effects. It is not difficult to imagine new sources 
of CP violation. In supersymmetry, for example, there are even too 
many. The problem is that we must find a model for new physics such 
as to obtain a sizeable contribution to e' / e while remaining within the 
stringent constraints imposed by e and by other measured quantities. 
This problem can be circumvented, so that for instance supersymme- 
try is potentially able (with some special assumptions) to produce the 
required effect on e'/e still fulfilling the other phenomenological con- 
straints |36| . 



At present, our preferred answer is the second one. Hopefully, improve- 
ments in non-perturbative techniques and a further insight in kaon phe- 
nomenology will clarify the mechanism responsible for the "large" value of 
e'/e and its connections with the AI = 1/2 rule. 
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